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I. INTRODUCTION 



Quantum mechanics gives predictions in form of probabilities. Already some of thc fathers of the theory were puzzled 
with question whether there can exist a deterministic structure beyond quantum mechanics which recovers quantum 
statistics as averages over "hidden variables" (HV) . In his famous impossibility proof Bell made precise assumptions 
about the form of possible underlying HV structure that allows mutually distant systems to be independent of one 
another [| . He derived the inequality which must be satisfied by all such (local realistic) structures and presented the 
example of quantum predictions which violate it. In this way the famous Einstein-Podolsky-Rosen (EPR) paradox 
was solved [2j. 

Noncommutativity of quantum theory precludes simultaneous deterministic predictions of measurement outcomes 
of complementary observables. For the EPR this indicated that "the wave function does not provide a complete 
description of physical reality" . They expected the complete theory to predict all possible measurement outcomes, 
prior to and independent of the measurement (realism), and not to allow "spooky action at a distance" (locality). 
Such a completion was disqualified by Bell. 

A more general version of two-particle Bell's theorem was given by Clauser, Horne, Shimony, and Holt (CHSH), and 
still extended by Clauser and Horne (CH) The important feature of the CHSH and CH inequalities, which hold 

for all local realistic theories, is that they can be not only compared with ideal quantum predictions, but also with 
experimental results. Thus a debaté that seemed quasi-philosophical, could be tested in the lab (see [5j for review of 
early experiments). 

The three or more particle, or as we now say qubit, versions of Bell's theorem were presented by Greenberger, Horne, 
and Zeilinger (GHZ), surprisingly 25 years after the original Bell's paper In contradistinction with the two particle 
case, now the contradiction between local realism and quantum mechanics could be shown for perfect correlations. 
Immediately after that Mermin produced series of inequalities for arbitrary many partides, which cover the GHZ 
case, and made the GHZ paradox directly testable in the laboratory Q. A complementary series of inequalities was 
introduced by Ardehali y|. In the next step Belinskii and Klyshko gave series of two settings inequalities, which 
contained the tight inequalities of Mermin and Ardehali Finally the full set of tight two setting per observer 
iV-party Bell inequalities for dichotomic observables was found independently in Refs. [ioL IÏÏL IÏ2T] . All these series of 
inequalities are a generalization of the CHSH ones 3] . Such inequalities involve only N party correlation functions. 

The process of finding new series of Bell inequalities (for many qubits, many settings per observer, involving 
all possible correlations - for arbitrary experimental arrangement) continues until this day, and is the tòpic of this 
overview. We shall try to present the current state-of-the-art, concentrating on multiqubit Bell inequalities for 
correlation functions. We will not give a detailed analysis of the assumptions behind Bell's inequalities. The reader 
can find them in excel·lent papers likeBEJ. 

With the emergence of the new sub-branch of physics (and information theory), Quantum Information, the Bell 
theorem, and Bell inequalities found applications far away from foundations of quantum physics. The security analysis 
of the first entanglement based quantum cryptography scheme involves Bell's inequalities 14]. This now is strengthen 
by the analysis of Scarani and Gisin, who showed that violation of Bell inequality is indecd a vàlid security criterion 
in quantum crypto-key distribution [Ï5|. Recently it was shown that with every Bell inequality one can associate 
a specific quantum communication complexity problem. With the use of quantum states which violate such Bell 
inequality one can always construct a quantum communication complexity problem which outperforms all possible 
classical ones This result, as well as the one of Scarani and Gisin, as well as many other ones, suggest that 

violation of Bell inequality is a criterion of direct usefulness of entanglement in quantum information processing |17| . 
The authors of this overview think that there are very many open qüestions associated with future generalisations of 
series of multiqubit Bell inequalities, and that it is still a fascinating field of studies, which will find new applications, 
and new surprising results. 
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II. BELL AND CLAUSER-HORNE-SHIMONY-HOLT 

The whole history of the EPR paradox the reader can find in the beautiful review by Clauser and Shimony [jj. 
Twenty nine years after the EPR paper Bell proved that the completion of quantum mechanics expected by EPR, 
is impossible [lj . In his original proof Bell utilized the perfect anticorrelations which arise whcnevcr Alice and Bob 
measure local spins (with respect to the same direction) on the two-qubit system in the state: 

^-) = ^[|0) A |1) B -|1)^|0) B ], (1) 

where |0) and |1) denotes the eigenbasis of the local <j z operator. However unavoidable experimental imperfections 
imply that correlations are never perfect. Here we re-derive the CHSH inequality, for which perfect correlations do not 
have to be assumed and thus the inequality can be directly experimentally checked. The violation of CHSH inequality 
implies that no local realistic explanation for the observed correlations is possible. But what if the inequality is 
satisficd? Can we then build a local realistic model for our observations? The answer is negative. Necessary and 
sufficient condition for local realistic model involves a set of inequalities, not a single one. 

The pair emission begins an experimental run. Alice and Bob measure one of two alternative settings each run. 
Their choices what to measure are absolutely free, uncorrelated (statistically independent) with the workings of the 
source. According to realism all possible measurement outcomes exist prior to and independent of the measurement. 
Moreover locality assumes that the outcomes of Alice depend on her setting only, and the same for Bob. We denote 
their predetermined local realistic results as A\, for Alice, and B\, B2 for Bob. For example if Alice chooses 
to measure setting "1" she obtains the outcome A\, if she chooses to measure "2" she obtains A 2 . Under realism 
assumption all possible measurement outcomes are defined, even if only some of them are actually measured. The 
experiments on qubits can give one of two results, to which we ascribe numbers, +1 and —1, i.e. A^,Bi = ±1. We 
form a "vector" out of the predetermined results of each observer: A = (A±, A2) and B = (B\, B2) in this case. One 
can also define a "vector" (or if you like, a "tensor") of local realistic correlation functions, Elr, with components 
E LR (k, l) = (AkBi) avg , where the average is taken over many experimental runs. Then all such local realistic models, 
Elr, can be written as: 

È L r= J2 P{A,B)A®B, (2) 

A,B=(±1,±1) 

where P{A,B) is the probability with which a certain quàdruple of predetermined results {A\, A2, B\, B2} appears. 
That is, every local realistic model of the correlation functions is a convex combination of the extreme points A® B, 
and thus lies within a convex polytope, spanned by the vertices A ® B. The necessary and sufhcient condition for 
local realistic description is a set of inequalities which define the interior of the polytope and are saturated at the 
border hyperplane of it. Such inequalities are called tight Bell's inequalities. 

Let us present a construction of the necessary and sufficient condition for the possibility of local realistic description 
of correlation functions for such experiments with two qubits. First we derive a necessary condition for local realistic 
model, then construct such a model proving that the condition is also sufficient. We introduce a more elaborated 
notation for future use in the generalisation to arbitrary number of partides. The two local dichotomic observables are 
parametrised by vectors n\ and ri 2 , and the predetermined results for the jth party by Aj{n\) = ±1 and Aj(n J 2 ) = ±1, 
as for now j — 1,2 (1 for Alice, 2 for Bob). Since Aj(n J k ) = ±1, for each observer j one has either |A,(n{) + 

^(^2)1 = an( i lAK^i) — Aj'Í^DI = 2, or vice versa. Therefore, for all sign choices of Si,S2 = ±1 the product 
[Ai(n}) + siAi(n\)][A2{n\) + S2A 2 (n 2 )] vanishes except for one sign choice, for which it is equal to ±4. If one adds 
up all such four produets, with an arbitrary sign in front of each of them, the sum is always equal to the value of the 
only non-vanishing term, i.e., it is ±4. Thus the following algebraic identity holds for the predetermined results: 

^I2,i2;s= J2 S(s 1 ,S2)[A 1 (nl) + s 1 A 1 (nl)}[A 2 (nl) + s 2 A 2 (n 2 2 )}=±4:, (3) 

Si ,S2=±1 

where S(sx, S2) stands for an arbitrary "sign" function of the summation indices s±, S2, such that its vàlues are only 
±1. The notation Ai2,i2;S describes the situation where two parties choose between two settings "1" or "2". 

After averaging the expression © over the ensemble of the runs of the experiment one obtains the following set of 
Bell inequalities: 

J2 S( Sl ,s 2 ) J2 ^i 1 ~ 1 4 2 ' 1 E(k 1 ,k 2 )\<A. (4) 

Sl,S2 = — 1,1 fei,fe2 = lj2 
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Since there are 16 different functions S(s\,S2), the inequalities Q represent a set of 16 Bell inequalities for the 
correlation functions. A specific choice of the sign function, 5(51,52) = v2sin(f'ïï' + (si + s 2 — 2)^), leads to the 
well-known CHSH inequality: Q 

\E(Í, 1) + E(í, 2) + E(2, 1) - E(2, 2)| < 2. (5) 
The set of all 16 inequalities l@J is equivalent to a single general Bell inequality: 

E I E ^ 1 - 1 4 2 " 1 ^i,fe)|<4. (6) 

si,s 2 =-l,l fcl,fc 2 = l,2 

The equivalence of JïïJ) and 10} is evident, once one recalis that, for real numbers \a + b\ < c and \a — ò| < c if and only 
if \a\ + \b\ < c, and writes down a generalization of this property to sequences of an arbitrary length. 

Whenever local realistic model exists inequality ijïïjl is satisfied by its predictions. To prové the sufhciency of 
condition © we construct a local realistic model for any set of correlation functions which satisfy it, i.e. we are 
interested in the local realistic models £7 (fe, fe) such that they fully agree the measured correlations E(ki, fe) for 
all possible observables fe, fe = 1,2. Recali that the set of local realistic correlation functions can be put as (2), 
with A — (Ai(n\), siAi(n\)) and B = (A 2 (n1), 52^2(712)). Let us ascribe for fixed Si,S2, a hidden probability that 
Aj(n{) = SjAj(n 3 2 )) (with j — 1, 2) in the form familiar from Eq. 

P{SUS2) = \\Y. S l~ ls 2~ lE ^ k z)\- ( 7 ) 

Obviously these probabilities are positive. However they sum up to identity only if inequality JHJ is saturated, otherwise 
there is a "probability dèficit", AP. This dèficit can be compensated without affecting correlation functions. 

First we construct the following structure, which is indeed the local realistic model of the set of correlation functions 
if the inequality is saturated: 

E £(si,s 2 )P(si,s 2 )(Mi)<8 (l,fla), (8) 

Si ,S 2 = ±1 

0. Now if AP > we add a "taü" to this 

(9) 

This "tail" does not contribute to the vàlues of the correlation functions, because it represents the fully random noise. 
Thus, the sum of © is a vàlid local realistic model for È = (E(l, 1), E(l, 2), E(2, 1), £7(2, 2)). 

Let us make some additional remarks. In the four dimcnsional real space where both Elr and E are defined one 
can find an orthonormal basis set S SlS2 = |(1, si) <8> (1, S2). Within these definitions the hidden probabilities acquire 
a simple form: 

P(sx,s 2 ) = ^\S SxS2 -È\, (10) 

where the dot denotes the scalar product in 1Z 4 . Now the local realistic models, Pla, can be expressed as: 

É LR = E \S S1S2 ■ ÈlA^nDA^nDSs^. (11) 

Si ,S 2 = ±1 

The modulus of any number \x\ can be split into |a;| = x sign(ï), and we can always demand the product Ai(n})j4 2 (rtJ) 
to have the same sign as the expression inside the modulus. Thus we have: 

È= E (Ss lS2 -É)S SlS2 . (12) 

si,s 2 =±l 

The expression in the bracket is the coefficient of tensor E in the basis S SlS2 . These coefficients are then summed 
over the same basis vectors, therefore the last equality appears. 

In this way the set of inequalities , or its equivalent - the single inequality I© ~ is proven to be sufficient and 
necessary for the possibility of local realistic description of correlation experiments on two qubits, where both Alice 
and Bob measure one of two local settings. This kind of reasoning can also be applied to arbitrary number of qubits. 



where S(si,S2) is the sign of the expression within the modulus in Eq. 
expression given by: 

AP 1* - 

"TeT £ 

A,B=(±1,±1) 
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III. MANY QUBITS 



Exiting features of the GHZ states led to a ràpid development of new Bell inequalities for multiqubit systems. Series 
of Bell inequalities for correlation functions were discovered by Mermin, Ardehali, Belinskii and Klyshko 0, S 0- 
Below we present a derivation of series which form the complete set of inequalities for TV parties, two-settings problem. 

The generalisation of the approach presented for two-qubit case to many qubits is straightforward. For N partides 
the generalisation of idcntity 10 consists of the sum of the products Aj(n J 1 ) + SjAj(n J 2 ) = ±2 for the jth party, and 
the summation is taken with more general "sign function", S(si, ...,sjv), of N parameters: 

N 

•lli ,•>:> S(s 1 ,...,s N )l[[A j (n{)+ Sj A j (4)]=±2 N , (13) 

si,...,sjv=±1 3 = 1 

Sincc thcre are 2 different sing functions S, the above formula leads to the set of 2 Bell inequalities. Using the 
trick described above we can write a single inequality equivalent to set of all 2 2 inequalities [Ï(Í IllL Il2j | : 

J2 I E s k ^-\..s k N -- 1 E(k 1 ,...,k N )\<2 N . (14) 

S 1 ,...,S N = — 1,1 fei,...,fcjy=l,2 

Many of these inequalities are trivial, e.g., if S(s\, s^) = 1 for all arguments, we get the condition \E(Í, 1, ...1)| < 1. 
Specific other choices give non-trivial inequalities. For example, for N = 3, we can recover the inequality (equivalent 
to the GHZ argument) given by Mermin: 7] 

|£(2,1,1) + £(1,2,1) + £(1,1,2)-£(2,2,2)| <2. (15) 

Up to now we have shown that if a local realistic model exists, the general Bell inequality (|14fl follows. The converse 
is also true: whenever inequality (|14fl holds one can construct a local realistic model for the correlation function, in 
the case of a Standard Bell experiment. For N partides the hidden probability that the predetermined outeomes of 
the jth observer are Aj(n{) — SjAj(n J 2 ) is given by the form familiar from Eq. l|14f) : 



P( Sl ,...,s N ) = ^ J2 s kí - 1 ...s k /- 1 E(k u ...,k N ) . (16) 



k\ . . . . , fcjv — 1,2 



The same steps as for two qubits above (now in the 1Z 2 space) lead to the result that any correlation experiment 
satisfying (|14|l can be explained within local realistic picture. This establishes the general Bell inequality (|14J) as a 
necessary and sufficient condition for local realistic description of N particle correlation functions in Standard Bell-type 
experiments. That is one can claim that the set of Bell inequalities represented by (|14fl is complete. This completeness 
implies that all series of Mermin A-qubit inequalities, which give tight inequalities, are a subset of the inequalities 
generated by l|14|). This also applies to tight Ardehali inequalities and the full set of Bclinskii-Klyshko inequalities 
E0. 



IV. MORE THAN TWO SETTINGS 



A general way to establish a necessary and sufficient condition for local realistic description is to define the facets 
of a correlation polytope 0, 120I |2Ï| . However this is a computationally hard NP problem • Even the complexity 
of printing out all Bell inequalities grows rapidly with the number of partides or settings involved! Now we present 
an efhcient method for generation of tight Bell's inequalities, which however do not form a complete set. This method 
was invented by Wu and Zong (2jJ |24| , and generalized in Ref . [2í| . 

We start with the case of N = 3 observers. Suppose that the first two observers can choose between four settings, 
and the third one between two settings. We denote such problem as 4 x 4 x 2. We already know that the local realistic 
vàlues satisfy the following algebraic identity: 

Ai2,i2,S'= E S '( Sí ' s *)[Miï\) + s 1 A l {nl)}[A 2 {nl) + S2A 2 {nl)} = ±4, (17) 

Sl ,S2=±1 

where S'{sx,S2) is any sign function, i.e. such that S'(sx,S2) = ±1. In analogous way, we can define ^34,34,5" by 
replacing A\(n\), A\{n\), A2(n\), A2{n\) by A\{n^), Ai(n\), A 2 (n 2 i ), ^(nf), respectively, and S' by S". Depending 
on the value of s = ±1 one has (Ai 2y s' + s^34,S") = ±8, or 0. By analogy to l|17l) one has: 



^1234,12= E S(si,s 2 )[A 12 , s > +s 1 A3^s"}[A 3 (n 3 1 ) + s 2 A 3 (nl)}= ±16. 

si ,s 2 =±l 



(18) 
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After averaging over many runs of the experiment, and introducing the correlation functions E(i,j,k) = 
(A 1 (n^)A 2 (n|)A 3 (íï|)) a vg one obtains multisetting Bell's incqualities. Because of the freedom to choose the sign 
functions S,S',S", we have (2 4 ) 3 = 2 12 Bell's inequalities. 

All these inequalities can be reduced to single "generating" inequality, in which all the sign functions S, S' , S" are 
non-factorable. The choice of factorable sign function is equivalent to having a non-factorable one, and some of the 
local measurement settings equal j2^] . The "generating" Bell's inequality can be chosen as (here all sign functions are 
equal to V2sin(§7r + {s 1 + s 2 - 2)f )): 



+ 
< 4 



A 3 (nl) + A 3 (n 3 2 )\ [^(^[(^(n?) + A 2 (n 2 2 )} + A x {fi\)[A 2 {nl) - A 2 {n\)] 
A z {n\) - A 3 (n 3 2 )] U (4)[A 2 (n§) + A 2 {n\)] + Ai(fíà)[A 2 (fí§) - A 2 (n\)\ 



(19) 



avg 



All other inequalities can be obtained by changes Aj {n° 



k j--Aj(fti) 

The method can be generalized to various choices of the number of parties and the measurement settings. Now we 



present the 2" 1 x 2 



)jV-l 



x 2 



N-2 



x 2 case. Consider N — 4 observers. We start with the identity l|18|) . One 



can introduce a similar formula for the settings {5,6, 7,8}, for the first two observers, and {3,4}, for the third one. 
The fourth observer chooses between two settings with local realistic vàlues A±{fi\) and A±(n\). Applying the same 
method as before, one obtains an identity which generates Bell's inequalities of the 8x8x4x2 type: 

^ S(si, s 2 )[Ai 234t i 2 + sxA^rs^A^nf) + s 2 A4(ft 2 )} =±64, 

Sl ,s 2 =±l 



(20) 



where ^1234,12 and j4s678,34 depend on some three sign functions. One may apply this method iteratively, increasing 
the number of observers by one, to obtain inequalities involving exponential (in N) number of measurement settings. 
As another example we construct the inequalities involving N partners, where first N — 1 observers choose one of 4 

settings and the last one chooses between 2 settings. We use the local realistic quantity Ai 2 12 dcfined in Eq. 

for N — 1 parties choosing between 2 settings each: 



112,...,12 



E 

,S N _i = ±l 



S'(s u ...,s N ^)l[[A j (nÍ)+s j A j (nÍ)}=^ 
3=1 



N-l 



(21) 



and analogically introduce ^34,... ,34 for another pair of observables available to each party. The sign function in 

A34 34 can be different from S'(si, sn-i) in Eq. (|21|l . By including the iVth observer, who can choose between 

2 measurement settings, we obtain: 



E 

Sl ,Si— ±1 



S(s!, fl 2 )(Ai2,...,ia + 8iAm,...,s*)(A n (H?) + s 2 A N (n 2 y )) = ±2 



loM+1 



(22) 



One can use this expression for generating Bell inequalities for N observers in the same way as it was previously done. 

In order to show the full strength of the method our next example gives a family of Bell inequalities for N — 5 
qubits, which involves eight settings for first two observers and four settings for the other three. We take the identity 
^4-1234,12 defined in (JTHJ, vàlid for the 4x4x2 case of three observers, and define a similar quantity for another set 
of 4 x 4 x 2 observables, namely ^5678,34- Note that the sign functions entering ^5678,34 can be different from those 
entering ^1234,12- For the other two observers we introduce: 



E 

si ,s 2 =±l 



5'(si, s 2 )(£>i + «iAOOEi + s 2 E 2 ) = ±4 



and a similar expression for another pair of observables D 3 , D 4 and £3, E4: 

^34,34 = 5 "( Sl > S2 )( J °3 + S 1^4)(£ 3 + S2-E4) = ±4- 



(23) 



(24) 



Sl ,s 2 =±l 



In the next step we get the following algebraic identity which can be used, via averaging, to generate a family of 2 28 
Bell inequalities: 



^ S(si, S 2 )(Ai 234) i 2 + SiA 5678>3 4)(yli 2i i 2 + S 2 ^34,34) = ±256. 



(25) 



si ,s 2 =±l 
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It is clear that there is no bound in extending this type of derivations. Finally let us note that all the inequalities 
with lower number of settings can be obtained from our construction by making some of the local settings identical. 

The multisetting inequalities constructed by the above procedure are tight. Consider the case of 4 x 4 x 2 in- 
equalities. The left hand side of the identity (|18|l is equal to ±16 for any combination of predetermined local 
realistic results. In a 32 dimcnsional real space, one can build a convex polytope, containing all possible local 
realistic models of the correlation functions for the specified settings, with vertices given by the tensor products of 
v = (Ai(ríí),Ai(n^),Ai(n^),Ai(n^)) ® (A 2 (n 2 ), A 2 (fí$), A 2 (nf), A 2 (nf)) <g> (A 3 (nf),A 3 (nD). It has 256 different ver- 
tices. Tight Bell inequalities define the half-spaces in which is the polytope, which contain a face of it in their border 
hyperplane. If 32 linearly independent vertices belong to a hyperplane, this hyperplane defines a tight inequality. Half 
of the vertices give in H18|) the value 16 and another half give —16. Every vèrtex, v from the first set has a partner 
—i) in the second one. Next notice that any set of 128 vertices v, which does not contain pairs v and —v contains 
a set of 32 linearly independent points (basis). Thus, each inequality is tight. This reasoning can be adapted to all 
inequalities discusscd hcrc. 

Most importantly, the multisetting inequalities reveal violation of local realism of classes of states, for which Standard 
inequalities, with two measurement settings per side, are satisfied. 



V. QUANTUM VIOLATIONS 

Bell inequalities are interesting only if there exists a certain state which violates them (for certain measurement 
settings). Let us derive the condition for violation of the general Bell inequality involving two measurement settings, 
1)14)1 . by an arbitrary two-qubit state. This result is a generalisation and reformulation of the condition given for two 
qubits by the Horodecki Family |2(j . We will use a decomposition of general mixed state of N qubits in terms of the 
identity operator ao = I and the Pauli operators Oi for three orthogonal directions i S {1, 2, 3}, given by: 

1 3 

P= 2 W T kx...k N <rki®-®<rk N , (26) 

/ci,...,fcív— 

where all the operators act in the Hilbert space of individual qubits. The (real) coefEcients T^...^, with kj = 1, 2, 3, 
form the correlation tensor T. Note further that T]^...^ = Tr(pak 1 <8 ■■■ &> &k N ), i-e. all coeficients are directly 
experimentally accessible. 

The full set of inequalities for the 2x2 problem is derivable from the identity l|17|) where we put non-factorable 
sign function: 



(A x + A 2 )B 1 + (A 1 - A 2 )B 2 ) <2. (27) 

/ avg 

All other inequalities are obtainable by all possible sign changes Xk — > — Xk (with k = 1,2 and X = A,B). The 
quantum correlation function E(dk,bi) is given by the sealar product of the correlation tensor T with the tensor 
product of the local measurement settings represented by unit vectors dk <S> &z, i.e. E(dk,bi) — (dk <8 bi) ■ T. Thus, 
the condition for a quantum state endowed with the correlation tensor T to satisfy the inequality (|77jl . is that for all 
directions a± , a 2 , b\ , b 2 one has 

( J ^- i ) ® ò i + Í- L· 2~ 1 ) ®b 2 \-T < 1, (28) 

where both sides of (|27|l were divided by 2. 

Note that the pairs of local vectors define the "local measurement planes" . Here we shall find the conditions for 
(1281) to hold for two, arbitrary but fixed, measurement planes, one for each observer. Therefore only those components 
of T are relevant which describe measurements in these two planes. Thus T is effectively described by a 2 x 2 matrix, 
or tensor T". 

Notice that À± = \{ax± d 2 ) satisfy the following relations: À + ■ À- = and ||A + || 2 + ||A_|| 2 = 1. Thus À + + À- 
is a unit vector, and A± represent its decomposition into two orthogonal vectors. Thus if one introduces unit vectors 
a± such that A± = a±a±, one has a\ + a 2 _ = 1. Thus one can put l|28|l into the following form: 

|S'·f'|<l, (29) 

where S = a + a + <g> b\ + a_a_ ® b 2 . Since a + ■ a_ = 0, one has S ■ S = 1, i.e. S is a tensor of unit norm. Any tensor of 
unit norm, U, has the following Schmidt decomposition U = \iï>i <8> Ü?i + X 2 v 2 ® w 2 , where Vi ■ Vj = 5ij, Wi ■ Wj = 
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and Af + A| = 1. The (complete) freedom of the choice of the measurement directions b\ and b 2 , allows one, by 
choosing b 2 orthogonal to b\, to put S in the form isomorphic with U. The freedom of choice of a\ and a 2 allows A+ 
and A_ to be arbitrary orthogonal unit vectors, and a+ and a_ to be also arbitrary. Thus S can be equal to any 
unit tensor. Therefore to get the maximum of the left hand side of <|28|) we put S = ^ 4^ f, and the maximum is 

||f'|| = Vf> ■ fi. Thus, 

max [ Y, T kl\ < 1. (30) 

fe,í=l,2 

where the maximization is taken over all local coordinate systems of two observers, is the necessary and sufficient 
condition for the inequality (|14|l to hold for quantum predictions. 

For 2 Ar ~ 1 x 2 JV ~ 1 x ... x 2 multisetting inequalities we can derive similar simple conditions for larger number of 
qubits. Consider the case of three qubits. For this situation we have the following inequality: 

|( Y / S(s u S2)[A 12iS ,+s 1 A Ui s,>}[C 1 + s 2 C2}) <16, (31) 
I \ ' I avg 

where S,S',S" are some non-factorable sign functions. The three-qubit quantum correlation functions E(Si,bj,Ck) 
can be represented as (oj ® Cj ® Cfc) ■ T (with the same meaning of the symbols as before; T is now a three index 
tensor). Thus the condition for the 4x4x2 inequalities to hold, in the quantum case, transforms into 

I [ÀX2,S> ® (ci + c 2 ) + À 3 4,s<< ® (ci - 82)] ■ f | < 8, (32) 

where e.g. 

À 12 , S ' = S'{s 1 ,S2)(a 1 + s 1 a2)<E){b 1 + s 2 b 2 ). (33) 

Sl,S2=±l 

To write down (|32|l we have used the freedom of introducing the sign changes Xi — > —Xi, compare (|19|) . By dcfining 
5(61 ±0*2) = c-(-c± which have the similar properties as a± and a±, inequality (|32(l transforms to: 

\c+Ài3,s> ®c+-f + c_i 3 4,s" ® c_ • f | < 4. (34) 
One can always choose c+ and c_ that maximize the left hand side. Since + c 2 _ — 1 this leads to the condition: 

[Ài2, s > ■ T (+) ] 2 + [À UiS » ■ f^] 2 < 4 2 , (35) 

where f^ is defined by = X)fc=i( c ±)fe^íifc> where in turn (c±)fc is the fc-th component of vector c±. Note that 
since c+ and c_ are orthogonal and normalized this procedure amounts to fixing of two (new) Cartesian axes for the 
third observer, and accordingly transforming the correlation tensor. Since À^s' depends on different vectors than 
^34, S" i one can maximize the two terms separately. Furthermore, since the problem of maximization of À nm g ■ 
is equivalent to the 2x2 case studied earlier, the overall maximization process gives the following necessary and 
sufficient condition for quantum correlations to satisfy the inequality (|31() : 



7 ' 2 

x=l,2 k-,l~=l,2 



max > ^ > j T 2 < 1. (36) 



When compared with the sufficient condition for 2x2x2 inequalities to hold, namely: 01 

max [ ]T T 2 lm ] < 1, (37) 

the new condition is more demanding because the Cartesian coordinate systems denoted by the indices ki,li and 
k 2 , l 2 do not have to be the same. Note further that if one uses in both terms in H36JI the same planes of observations 
(for the first two observers) then the condition reduces to (|37|l . which now is a sufficient and necessary one (but for 
4x4x2 settings). Thus the condition is necessary and sufficient whenever the vectors dcfining measurements for 
each party are limited to one plane. 
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TABLE I: The examples of necessary and sufficient conditions for violation of multisetting inequalities. 



N 


case 


Cn (the condition) 


2 


2x2 


Z)fc,!=l,2 T kl - 1 


3 


4x4x2 




=1,2 T fei2 + Z)v,!'=l,2 T fc'!'l - 1 


4 


8x8x4x2 




=1,2 ^*ili22 + Sfc 2 ,i 2 = l,2 71 fe 2 Í2l2 + 








1,2 T k 3 l 3 2\ + Sfc 4 , i4 =l, 2 T fe 4 Í4ll ^ 1 



In a similar way one can reach analogous conditions for violation of 2 N ~ 1 x 2 JV ~ 1 x 2 2 x ... x 2 inequalities by 
quantum predictions. In the Table 1 we present these conditions for small N. One can note a useful recurrence that 
can be used to write down the condition for arbitrary N. Let us define: 

E T li- ( 38 ) 

k,l=l,2 

Then the condition for two qubits reads: max(C2) < 1. Next let us put a recursive defmition: 

Cn = [Cjv-i]©2 + [Cjv-i]©i, (39) 

where [C^v— is the expression in the condition for N — 1 qubits in which the correlation tensor elements Tj 1 ...j íf _ 1 
are replaced by Tj 1 ...j JV _ 1 fc, i.e. elements of the AT-qubit correlation tensor. The "prime" denotes the fact that the 
second term does not involve components of T in the same set of coordinate systems (for the first N — 1 observers) 
as the unprimed term. 

The sufficient and necessary condition for N qubits to satisfy all 2 N ~ 1 x 2 N ~ 1 x 2 N ~ 2 x ... x 2 inequalities, within 
this convention reads: 

max(CAr) < 1. (40) 



VI. GISIN'S PROBLEM 



The theorem of Gisin states that any pure non-product state violates local realism, i.e. there are sets of measure- 
ments that can be performed on the state which cannot be described within local realistic picture [27l l2Sj . This 
theorem formalizes the intuition that entanglement is a purely quantum phenomenon. Using the approach presented 
here we can write down the following proof of Gisin's theorem for two qubits. Any state of two qubits is given in its 
Schmidt basis by — cosa|00) + sina|ll), with a £ [0, 7r/4]. The correlation tensor of this state has the following 
coefficients: T xx = sin2a, T yy = — sin2a, T zz = 1. Therefore the necessary and sufficient condition for local realism 
is violated for all non-product (a ^ 0) states: 

E ^ = l + sin 2 2a>l. (41) 

k 1 l—{x,z} 

As we shall show here the series of the two-settings inequalities (|14|l for TV qubits are failing to show violation of 
local realism for an important class of entangled pure states. This looks, superficially, as a counterexample to Gisin's 
theorem. But the theorem uses all possible measurement scenarios, which is not the case in the case of inequalities l|14|l . 
Furthcr, as it will be shown, the more-than-two setting inequalities do show violations for this class of states. This is 
a strong argument for continuation of the research towards finding inequalities with even more general structure. 

Scarani and Gisin noticed a surprising feature of the following state: 29] 

|V>) = cosa|0...0) + sina|l...l>. (42) 

They showed that for sin 2a < 1 /VW^ the states do not violate the Mermin-Ardehali-Belinskii-Klyshko (MABK) 
inequalities 0, @- This has been numerically obtained for = 3,4, 5 and conjectured for A^ > 5. Their result 
contrasts the case of two qubits and is highly counterintuitive as the states l|42f) are generalization of the GHZ states, 
which violate maximally the MABK inequalities |(| . The problem of Scarani and Gisin was later analysed using the 
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general two-qubit correlation function Bell inequality for arbitrary number of qubits l|14fl 30] . The family (|42|) does 
not violate any such Bell inequality if a and N are such that sin2a < l/\/2 N ~ 1 and N is odd. However for even 
number of qubits the general Bell inequality imposes stronger constraints for local realistic description, than MABK 
inequalities. 

What are the reasons for the completely different behaviour for N even and N odd? The expression 
Y)j (1 k N =xyTk 1 fcjvi which appears in the condition for violation of two-setting inequalities, can be understood 
as a "total measure of the strength of correlations" in mutually complementary sets of local measurements (as de- 
fined by the summation over x and y) |3l| . Then the unity on the right-hand side of the condition is the classical 
limit for the amount of correlations. Specifically, pure product states cannot exceed the limit of 1, as they can show 
perfect correlations in one set of local measurement directions only. In contrast, entangled states can show perfect 
correlations for more than one such set [IJ. Now, only if N is even, the states l|42[) already show perfect correlation 
between measurements along z-directions (as the product is then always +1) reaching therefore the classical limit 
(we assume a z eigenstates give the qubit computational basis). Yet, they also show additional correlations in other, 
complementary, directions. In the case of N odd, however, there are no perfect correlations along z-direction and the 
correlations in the complementary directions do not suffice to violate the bound of 1. 

Since the condition for multipartite two-setting inequalities to hold: 

max J2 T l-k N < 1, (43) 

fei,...,fcjv=l,2 

is only necessary, if a state satisfies it one cannot draw any conclusions about its possible local realistic model. 
Thus it could be that the states Ij42(l can violate local realism if we use some more optimal tool to test it. Indecd 
multisetting inequalities described extend the class of quantum states which do not admit local realistic explanation. 
Also the generalized GHZ states cannot be classically explained. Their nonvanishing correlation tensor components 
are (directions 1,2,3 are denoted by x,y,z; the basis {|0), |1)} is the eigenbasis of <r z ): T z ,.. z — cos 2a, for N odd, 
and 1 for N even, T x ... x — sin2a, and the components with 2k indices equal to y and the rest equal to x take the 
value (— l) k sin2a (there are 2 N ~ 2 such components). Let us assume that the last observer can choose only between 
settings x and z. Thus, we obtain for the condition for violation of multisetting Bell's inequality for odd number of 
observers (generalization of the condition (J2HI) 

V T 2 4- V T 2 

Z_-*fei ,...,fcjv_ i— x,y ki . . ix r k\ ,&ív- i—x,z fci.../cjv — \z 

= 2 W - 2 sin 2 2a + cos 2 2a > 1. (44) 

Thus, the Bell's inequalities are violated for the whole range of 7r/4 > a > and for arbitrary N in contrast to the 
case of Standard Bell's inequalities. 



VII. PROSPECTS 



We have described some series of Bell's inequalities for arbitrary number of qubits involving iV-particle correlation 
functions. In the case of two settings per observer the inequalities form a complete set, which gives the necessary 
and sufRcient condition for local realistic description of correlation experiments. Even for this simple scenario it is 
still an open question how would the necessary and sufRcient condition look like for more general case involving all 
correlations of all ranks between observers, such as those described by correlation functions for N — 1 parties, etc. 
(generalisations of CH inequality). The other problem is to consider more settings per observer. A step towards this 
direction was recently made by one of us, who constructed a full set of inequalities for three settings per observer and 
arbitrary number of parties 32]. The final aim is to construct such inequalities for arbitrary experimental setup. 

It is clear that only entangled states can violate Bell's inequalities, as any separable state has a local realistic model. 
For pure entangled states there always exists Bell's inequality that is violated by certain measurements performed on 
the state. But which mixed states violate local realism? Perhaps a step forward could be achieved by introduction of 
other, additional to Bell's (local realism) assumptions, based on the fundamental symmetries of physical laws, which 
would narrow the class of physically admissible local realistic theories [33, 01 • 

The work is part of the MNil Grant no. 1 P03 04927. MZ acknowledges Professorial Subsidy of FNP. WL and TP 
acknowledges a FNP stipend. 
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